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I. INTRODUCTION 



This thesis investigates the use of filtering techniques 
as an electronic counter countermeasure (ECCM) to overcome the 
effects of the optimum jammer operating against coherent phase 
shift key (PSK) and frequency shift key (FSK) receivers. These 
receivers are well-known solutions to the problem of optimum 
detection of known signals in the presence of additive white 
Gaussian noise (WGN) [Ref. 1] . The structure of such receivers 
is shown in Figure A.l. These two receivers represent the 
optimum signal processing algorithms for the discrimination 
of the digital signals s^(t) and Sq ( t) representing digital 
"one" and "zero" logical states respectively. To reduce the 
effectiveness of these optimum receivers. Ref. 2 and Ref. 3 
developed optimum jammers to be used against the PSK and FSK 
coherent receivers. From the resultant jammer waveforms derived 
in these references, this thesis investigates the use of front- 
end filters inserted in these receivers as shown in Figure A. 2 
as a method for reducing jammer effectiveness. 

From Ref. 2, the optimum power constrained jammer for 
either the coherent PSK or FSK receiver is a signal propor- 
tional to the correlator signal, shown in Figure A. 2 as s^(t). 
This correlator signal (see Ref. 1) is proportional to the 
difference of s-^(t) and s Q (t) , the digital "one" and "zero" 
signals. Though the correlator signal can take different forms 
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in these receivers depending on the modulation used, the theory 
used to develop these receivers shows that a jammer waveform, 
proportional to the difference of s^t) and s Q (t), is optimal. 

The jammer, therefore, takes the form of a continuous wave 
(cw) or tone centered at the frequencies of s^t) and s Q (t). 

In the case of the PSK coherent receiver, the jammer is a single 
tone centered at the signaling frequency. For the FSK case, 
the jammer waveform consists of the difference of two tones 
at the frequencies of s^(t) and Sg(t). The use of tones or 
cw jammers in angle modulated analog or digital receivers is 
discussed in Ref. 4. 

Given these forms of jammer signals, an obvious choice for 
the front-end filter is a band reject or bandstop filter. The 
filters chosen for analysis in this thesis are the ideal band- 
stop filter and a single zero second order bandstop filter. 
Though the ideal filter is not realizable, the analysis involv- 
ing its use will provide both insight into the problem and 
approximate results before investigating the use of the more 
complex second order filter. For the PSK coherent receiver, 
the bandstop regions of both the ideal and second order filters 
are centered at the frequency of s^(t) which is also equal to 
the frequency of Sg(t) . Since in this thesis a filter with 
only a single bandstop region will be analyzed, the bandstop 
region of the ideal and second order filters used in the FSK 
coherent receiver will be centered midway between the frequen- 
cies of s^(t) and Sg(t). 
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The analysis of the two forms will begin with the statis- 
tics of G, the output of the coherent receiver. From Ref. 2, 

G is a conditional Gaussian random variable with conditional 
statistics depending on whether s^(t) or Sq ( t) was transmitted. 

The analysis will then focus on determining the statistics of 
the random variable G', the output of the receiver with a front- 
end filter as shown in Figure A. 2. The mean and variance, 
conditioned on whether s-^(t) or Sg(t) was transmitted will 
then be determined using frequency domain analysis techniques. 
Finally, the general analysis will conclude with the determina- 
tion of an expression for the probability of error of the coherent 
digital receiver using a front-end filter to counter the opti- 
mum jammer. Next, the specific frequency domain forms of the 
PSK and FSK signals and jammers will be developed and used to 
calculate the probability of error performance for these re- 
ceivers with and without front-end filters. 

To determine the Fourier transforms of the jammers, this 
thesis will first assume that the jammer waveform is a time 
truncated function. Since in References 1, 2, and 3, the 
correlator signal s^(t) is defined over the interval (0,T), 
and the optimum jammer (from Ref. 3) is proportional to s^(t) , 
the optimum jammer will be treated as defined in the interval 
(0,T) also. Since from a practical standpoint, the jammer is 
a cw signal, the analysis will be modified in the sequel to 
account for this interpretation of the jammer. 

In Chapter II, a development of the mathematical expressions 
in the form of frequency domain equations for the statistics 
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of G' (the output of the coherent receiver of Figure A. 2) 
is undertaken. From these statistics, the expression for 
the receiver probability of error is developed. Chapter III 
demonstrates the application of these general expressions to 
the binary PSK (BPSK) coherent receiver problem with both 
ideal and second order front-end filters. Chapter IV presents 
similar analysis to that in Chapter III, applied now to the 
binary FSK (BFSK) coherent receiver. Chapter V presents 
quantitative and graphical analysis of the four receiver/ 
filter combinations, under the assumptions of both the cw and 
time truncated jammer forms. 



13 



II. ANALYSIS OF DIGITAL COHERENT RECEIVERS 
WITH FRONT END FILTERS 



In order to investigate the effects of filtering at the 
front-end of digital coherent receivers, it is necessary to 
review the operation and performance of digital coherent 
receivers in the presence of noise and jamming without the 
use of filters to notch out the jammer. As previously pre- 
sented in Ref. 2, the structure of the receiver to be analyzed 
is shown in Figure A.l. In the absence of a jamming signal 
this receiver is optimum for determining whether a "one" 
signal or a "zero" signal was transmitted in a given interval 
(0,T), with minimum probability of error [Ref. 1]. 

Following the analysis of Ref. 3, which assumes that the 
jammer is a deterministic waveform unknown to the receiver, 
the input signal r(t) to the receiver front-end is mathematically 
modeled as 

r (t) = s i ( t ) + n ( t ) + n^(t) 0 < t < T, i = 0,1 (2.1) 

where the s^(t) (i = 0,1) are used to transmit the digital "one" 
and "zero" data, n(t) is a sample function of a white Gaussian 
noise process with power spectral density level of N^/2 watts/ 
hertz, and as stated previously n^ (t) is the jammer waveform. 

From the analysis of Ref. 2, the output of the receiver is 
given by 
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G 



( 2 . 2 ) 



(r ' s d } 






where 



T 9 

(x,y) = / x(t)y(t)dt and ||x|| = (x,x) 

a 



Since G is a (conditional) Gaussian random variable, the sta- 
tistics are completely determined when the conditional mean 
and variance are found. The conditional mean is given by 

E [G | s^transmitted] = (s^,s^) + (n^,s^) 

+ i[||s o l| 2 -||s 1 l| 2 J (2.3) 

i = 0,1 

and the conditional variance is given by 

N 0 2 

VAR[G | (t) transmitted] = | | s^ | | , i = 0,1 . (2.4) 

A. FILTERING EFFECTS ON THE MEAN AND VARIANCE 

The effect of placing a filter at the front-end of the 
receiver is now analyzed by evaluating the conditional mean 
and variance of the receiver output. Figure A. 2 shows the 
placement of the filter at the receiver front-end. It can be 
seen that G becomes G’ where 
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C = (r',s d ) + |f| I s o I I 2 - ||s 1 || 2 ] 

(2.5) 

oo 

r' = r *h = / h(t-a)r(a)da 



and h(t) is the front-end filter impulse response. Since 

r ' (t ) = s^ ( t) + n'(t) + n!(t) i = 0,1 (2.6) 

it is simple to see that 

E [G' | s i (t) transmitted] = (s|,s^) + (n^s^) 

+ | |s o | | 2 - | | Sl | | 2 ) (2.7) 

i = 0,1 



where 



s|(t) = Si (t) *h(t) 

n ^ ( t ) = nj (t) *h(t) 

It is much more convenient to specify filters in the fre- 
quency domain because convolutions in the time domain become 
multiplications in the frequency domain. Thus looking 
specifically at the term (s|,s d ) where 

T 

(s!,s,) = / s!(t)s,(t)dt i = 0,1 (2.8) 

i a ' i a 
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or 



(s i' s d ) 



00 

/ s|(t) s dp ( t) dt i = 0,1 



(2.9) 



where 



s dp (t) = s d (t) ‘ p(t) 



with 



p(t) 



1, 0 £ t £ T 

0 , elsewhere 



and replacing s^^Ct) by its inverse Fourier transform equiva- 
lent (s d p(t) <=> S d p(oo) ) , namely 



s , (t) 
dp 



1 ( <-• , > j co t . 

— J S (m)e J dm 

2 77 _J dp 



( 2 . 10 ) 



yields 



(s!,s d ) 



OO 00 

/ s!(t)[^r / S (oj) e^ (ot dw] dt (2.11) 

l 2 tt __ 1 dp 



Interchanging the order of integration yields 
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(s i 



s d> - 



2tt 



OO CO 

/ Sdp (a,) 1 / si(t)e 3ut dt]dw i = 0,1 (2.12) 



and recognizing the fact that 



/ s!(t)e ja)t dt = S|(-a)) 



i = 0,1 



results in 




S d ) 




Sdp (w) (-0)) dcj 



i = 0,1 



(2.13) 



Furthermore, since 



s|(t) = s L (t) * h ( t) i = 0,1 

S^(o)) is given by 

S|((d) = S i (a))H(oj) i = 0,1 



and therefore 



<s i’ s d ) 



oo 

2~ / s dp (cj) S i (-w) H (-cj) dcj i = 0,1 



dp 



(2.14) 
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Replacing s^(t) with s dp (t) 



(n ! , sj becomes 
j d 



as in Equation 2.9, the term 



oo 

(n j ,S d ) = J n j (t) S dp (t)dt 



(2.15) 



and using Equation 2.10 in Equation 2.15 yields 




/ n'(t)[ 



2t t 



/ s (w) e jwt dw]dt 



(2.16) 



As done in Equations 2.11-2.13, interchanging the order of 
integration and using 



00 

N ! (-co) = / n! (t)e Da)t dt 

3 — OO 1 



(2.17) 



with 



NjU) = N_. (oj) H (tu) 



finally yields 




OO 

A / s, (U))N (-d))H(-a>)da) 

2tt ’ dp l 

— oo ^ 



(2.18) 
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The conditional mean of g ' , given that s^(t) is transmitted, 
is 



E[G'-|s^(t) transmitted] 



oo 

= 2¥ J S dp (m)s i (-u))H(-a J ) 



dto 



w 

+ Jr / S , (w)N . (-w) H (-to) da) 

2 it J dp j 



+ I l s 0 l I 2 ~ I |s x | | 2 ] 



(2.19) 



i = 0,1 



The variance of g' given that s^(t) is transmitted can be 
shown to take on the form 



VAR[G'|s^ transmitted] = E [ ( n ' ( t) , s^ ( t ) ) ] 



( 2 . 20 ) 



The right hand side of this equation becomes 



E [ (n 1 ( t) ,s d (t) PI 



= E [ / n'(t)s d (t)dt / n'(i)s d (x)di] 

— nr> — nn 



( 2 . 21 ) 



Equation 2.21 can now be written as 



E [ ( n ' (t) ,s d (t)) 2 ] 



OO 



/ 



— oo 



/ E[n' (t) n' (x) ] s d (t)s d (x)dtdx 



( 2 . 22 ) 
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Recognizing that the expected value operation in the integral 
is the autocorrelation function of the filtered white Gaussian 
noise, R n , (t-t) , Equation 2.21 can be written as 



E[(n',s d ) ] 



OO CO 



/ / V (t ‘ T)s dp <t)s d p (T)dtdT 



— oo — oo 



(2.23) 



Again, it is much more convenient to continue the analysis 
in the frequency domain. Using the inverse Fourier transform 
for the frequency domain representation of R n ,(t-x) and sub- 
stituting into Equation 2.23 yields 



E [ (n' ,s d ) ] 



OO OO 



/ _/ [37 / s n- <<*.)ei“ (t - T) d lu ] Sdp (t) 3 dp (T) 



—OO —OO 



dtdx 



(2.24) 



where s d ^(t) is defined by Equation 2.9. Interchanging the 
order of integration. Equation 2.24 becomes 



E[ (n 1 ,s d ) ] 



= 27 / s n ' (a))da) / s dp (t)eDa3tdt / s d P (T)e I|WTdT 



(2.25) 



It is easily recognized that 
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CO 



/ 



s , ( t) w ^dt 
dp 



S dp ( “ u) 



S* tw) 

dp 



and 



J s dp (T)e DajTdT = S dp (u)) 



Therefore in the frequency domain, VAR(G'|s^ transmitted) 
becomes 



oo 

VAR[G ' I s • ( transmitted] = / S , ( oj ) I S ( to ) I 2 du) (2.26) 

1 2tt 1 n 1 dp ' 

— oo *- 



Finally since n'(t) = n(t) * h(t), it is well known that 
S n' (a)) = s n (w>l H (w)| 2 

Since n(t) is a sample function of a random process that was 
assumed to have power spectral density S^w) = N^/2, for all cj , 



Var[G'|s^(t) transmitted] 



N °° 

^ / |H(u>) | 2 |S dp (ui) | 2 da) (2.27) 



2 2 

where the fact that | H (to) | |S (co) | is an even symmetric 
function has been used to reduce the limits of integration 
to half the real line. 
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In summary , G' is a (conditional) Gaussian random variable 



with conditional mean and variance given by 



E[G'|s^(t) transmitted] 



_ 1 _ 

2tt 



oo 




S dp ^ S j_ ( - w) H (- 00 ) doo 



JL f 

2 it J 



S dp (co)N j (-a))H(-w)dt0 + I [ 
i = 0,1 



( 2 . 28 ) 



and 



VAR[G’|s^(t) transmitted] 



N, 



= — ( 
2 7T -J 



H (co) 



■ s dp ( “> 



"doo 



i = 0,1 



( 2 . 29 ) 



B. PROBABILITY OF ERROR PERFORMANCE 

The analysis thus far has been concerned with determining 
the statistics of the random variable G', the output of the 
coherent receiver. The next step is to determine the proba- 
bility of error from the decision process, where the amplitude 
of G' is compared with a threshold y. Referring to Figure A. 2, 
if the value of G' is greater than y, the receiver decides that 
a "one" or s^ (t) was transmitted, and conversely if the value 
of G' is less than y the receiver decides that a "zero" or 
Sq ( t) was transmitted. The quantity, y, is determined in the 
derivation of the optimum receiver for detecting known signals 
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in additive white noise [see Ref. 1] . That analysis determines 
the threshold of the correlator (coherent) receiver as 



Y 




0 



(2.30) 



However the receiver of Figure A. 2 treats the second term of 
Equation 2.30 as a bias, thus the threshold becomes 



where Nq/ 2 is the power spectral density level of the white 
Gaussian noise process. Here, n is determined by the specific 
decision rule applied. In this case, the receiver strategy 
is to detect s^ (t) or Sg(t) with minimum probability of error, 
so from Ref . 1 



All further analysis in this thesis will assume that the 
probability that s-^(t) was transmitted is equal to that of 
Sq ( t) being transmitted, so the decision threshold is 



Now, the receiver can make two types of errors. Specif- 
ically, deciding that s Q (t) was transmitted when in fact s 1 (t) 



Y 




(2.31) 



P{s^(t) was transmitted} 



(2.32) 



n 



P { Sg ( t ) was transmitted} 
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Y 



0 



(2.33) 
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was transmitted (defined as error 1), and deciding that s^(t) 
was transmitted when s^Ct) was actually transmitted (defined 
as error 2). The receiver decides that Sq ( t) was transmitted 
when G' is less than the threshold y, and that s^(t) was trans- 
mitted when G' is greater than y . So, the probability of error 
1 occurring is the probability that G' is less than y when s^(t) 
is transmitted. Since G' is a Gaussian random variable with 
conditional probability desnity functions dependent on whether 
s^(t) or SQ(t) was transmitted 



Y 



(g'-ir^) 

2v 



error 1 



/2ttv. 



dg’ 



(2.34) 



where m^ = E[G'|s^(t) was transmitted] and v^ = VAR[G'|s^(t) 
was transmitted]. By letting x = g'-m^//v , Equation 2.34 
becomes the integral of the standard Gaussian probability 
density function also known as the error function (erf (x) ) 
of the form 
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error 
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-V /V g 
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/2tt 



-x 2 /2 



dx 



(2.35) 



The probability of error 2 occurring is the probability that 
G' is greater than y when s 0 (t) is transmitted. Therefore 



- (g ' -m Q ) /2v 



error 2 
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(2.36) 



where 



nig = E[G'|s^(t) transmitted]. 

Using a similar change of variables as in Equation 2.35 
yields 
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error 
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e dx 



(2.37) 



where y is zero from Equation 2.33. Finally, the average 
probability of error becomes 



P 

e 



P P error 1 + 



(1-p) P 



error 2 



(2.38) 



where p = P{s^(t) was transmitted} and 1-p = P{sQ(t) was 
transmitted}. Since p = 1/2 and substituting Equations 2.35 
and 2.37 into 2.38, the average probability of error is 
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(2.39) 



III. PHASE SHIFT KEY COHERENT RECEIVER 
WITH A FRONT-END FILTER 



A. PHASE SHIFT KEY WITH AN IDEAL FRONT-END FILTER 

Using the equations for the mean and variance of the con- 
ditional Gaussian random variable G' (Equations 2.28 and 2.29), 
along with the probability of error equation (Equation 2.39), 
it is possible to investigate the effect of a front-end filter 
on the probability of error performance of the coherent receiver. 
The receiver analyzed in this case is the optimum receiver for 
binary phase shift key modulated signals with an ideal band- 
stop filter in the front-end. The filter will be characterized 
as having a gain of unity over all frequencies, with the excep- 
tion of the regions from w = ±Wq -B/2 to w = ±u)q + B/2 where 
the gain is zero. The bandwidth of the filter is B radians/ 
second and the center frequency is u)q radians/second. Although 
this filter is not realizable, it serves both as a good approxi- 
mation to a higher order bandstop filter, and as a simple 
example so as to gain insight to the receiver performance. 

Using Figure A. 2 as a reference, the generic coherent 
digital receiver can be made into a BPSK receiver by setting 

s,(t) = 2A cos o) n t 0 < t < T (3.1) 

d U — — 

and modeling the input PSK signals as 
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s^t) = A cos(u> Q t + e^(t)) 0 <_ t <_ T, i = 0,1 (3.2) 

where 

! tt for i = 0 
0 for i = 1 

Thus s^(t) (i = 0,1) can be expressed as 
S Q ( t) = -A COS (Ji^t 

0 £ t £ T (3.3) 

s 1 ( t) = A COS U) Q t 

The jammer waveform is set to 

n j(t) = /P~ /27T cos w Q t 0 £ t £ T (3.4) 

since in Ref. 2 it was shown to be optimum against BPSK trans- 
mission. Here P n ^ is the jammer output power. Finally the 
noise input into the receiver of Figure A. 2 is a white 
Gaussian random process with power spectral density level of 
Nq/ 2 watts/hertz. 

1 . Calculation of the Variance 

From Equation 2.27, the variance of G' conditioned on 
s^(t) being transmitted is given by 
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VAR[G'|s^(t) transmitted] 

N °° 

= 27 / |H(0)) | 2 |S dp (o)) | 2 da) i = 0,1 



(3.5) 



Since the front-end filter is ideal. Equation 3.5 can be 
modified to become 



VAR[G'|s d (t) transmitted] 
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= — [ ( 
2 it l J 



oJq+B/2 



S dp (w) I da> “ 



a) 0 -B/2 



S dp (a)) I da) 



(3.6) 



i = 0,1 



since | H ( uj ) | is unity except over the bandstop region. 



In order to calculate |s^ (w)| , since s dp (t) = s d (t)p(t) 



where 



p (t) 



! 1 0 < t £ T 

0 elsewhere 



then 



s dp (t) = 2A p(t) cos p^t (3.7) 

The Fourier transform identity f(t)*g(t) <^> (1/2 tt) xf(oj) * G(to) 

will be used in order to calculate S, (to) . Since P(co) is 

dp 
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P(w) 



T Sa (toT/2) e 



(3.8) 



- jcoT/2 



and 



S d (w) = 2A7r [6 ((jj-o)q ) + 6(a)+a)g)] 

S, (oj ) is the convolution of these Fourier transforms, result- 
dp 

ing in 

-j (u-ojg) T/2 

S dp(“) = AT [Sa (03-w Q ) je 

-j (u+o) 0 ) T/2 

+ Sa (oj +co g ) 2& ] (3.9) 

2 

From this, | ( oj ) | is obtained by multiplying S^fw) by 

S, (-oj) which results in 
dp 

l s dp (w)| 2 = A 2 T 2 • [Sa^u-Wg)^ +Sa 2 (a)+Wg)|] (3.10) 

where the product of Sa (w+Wg ) T/2 and Sa (w-oig ) T/2 has been 
ignored because it is essentially zero for practical values 
of cog and T. 

Since the variance consists of an integral from oj = 0 
to infinity, for large ojg the contribution of the Sa(w+Wg)T/2 
term is negligible. Therefore, substituting into Equation 3.6 
yields 
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Var[G'|s^(t) transmitted] 



2 2 N n 00 2 T 

= A JL[ I Sa (oo~(jl) 0 ) idoo 



Wq+B/2 



Wq-B/2 



Sa^a)-a) 0 )|du)J (3.11) 



i = 0,1 



In this form Equation 2.11 is unmanageable due to its complexity 
and the number of variables present. To simplify, the follow- 
ing substitutions will be made: 



1. x = (oj-Wq) T/2 , 2/Tdx = d w 

2 . for oo = 0 , x = - 00 qT /2 which for large values of oo^T 
becomes approximately -oo 

3 . for oo = oo , x = oo 

4. for oo = oo q ± B/2 , x = +BT/4 

Equation 3.11 now becomes 



VAR[G'|s^(t) transmitted] 



A 2 TN n oo BT/4 2 

[ j Sa^x)dx - j Sa(x)dx] (3.12) 

11 -oo -BT/4 

i = 0,1 



2 

Finally the Sa (x) is a symmetric function about x = 0 and 

oo 2 

f Sa (x)dx = tt/ 2 from math tables. Thus the variance of G' 

0 
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